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Abstract 

In this paper, we continue to investigate the thermodynamic properties of stellar 
self-gravitating system arising from the Tsallis generalized entropy. In particular, 
physical interpretation of the thermodynamic instability, as has been revealed by 
previous paper (Taruya & Sakagami, Physica A 307 (2002) 185), is discussed in detail 
based on the framework of non-extensive thermostatistics. Examining the Clausius 
relation in a quasi-static experiment, we obtain the standard result of thermody- 
namic relation that the physical temperature of the equilibrium non-extensive sys- 
tem is identified with the inverse of the Lagrange multiplier, T p h ys = 1/(3. Using 
this relation, the specific heat of total system is computed, and confirm the com- 
mon feature of self-gravitating system that the presence of negative specific heat 
leads to the thermodynamic instability. In addition to the gravothermal instability 
discovered previously, the specific heat shows the curious divergent behavior at the 
polytrope index n > 3, suggesting another type of thermodynamic instability in the 
case of the system surrounded by the thermal bath. Evaluating the second variation 
of free energy, we check the condition for onset of this instability and find that 
the zero-eigenvalue problem of the second variation of free energy exactly recovers 
the marginal stability condition indicated from the specific heat. Thus, the stellar 
polytropic system is consistently characterized by the non-extensive thermostatis- 
tics as a plausible thermal equilibrium state. We also clarify the non-trivial scaling 
behavior appeared in specific heat and address the origin of non-extensive nature 
in stellar polytrope. 
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1 Introduction 



Due to its complexity and peculiarity, stellar self-gravitating system has long 
attracted much attention in the subject of astronomy and astrophysics, and 
even statistical physics. For an isolated stellar system, the dynamical equi- 
librium is rapidly attained after a few crossing time and the thermodynamic 
description provides useful information in characterizing the late-time behav- 
ior of this system. Even in this simplest situation, however, the equilibrium 
state of self-gravitating system shows various interesting phenomena, which 
may offer an opportunity to recast the framework of the thermodynamics 
and/or statistical mechanics. 

In earlier paper, applying the Tsallis' generalized entropy [1], we have stud- 
ied the thermodynamic instability of self-gravitating systems[2]. The self- 
gravitating stellar system confined in a spherical cavity of radius, r e , exhibits 
an instability, so-called gravothermal catastrophe, which has been widely ac- 
cepted as a fundamental physical process and plays an important role for the 
long-term evolution of globular clusters [3-5] . The presence of this instability 
has been long known since the pioneer work by Antonov[6] and Lynden-Bell 
& Wood[7]. Historically, the gravothermal catastrophe has been studied on 
the basis of the maximum entropy principle for the phase-space distribution 
function, with a particular attention to the Boltzmann-Gibbs entropy [8,9]. 

In contrast to previous work, we have applied the Tsallis-type generalized 
entropy to seek the equilibrium criteria for the first time (for comprehensive 
review of Tsallis formalism and its application to other field of physics, see 
[10,11]). Then, the distribution function of Vlassov-Poisson system can be re- 
duced to a stellar polytropic system[12,13]. Evaluating the second variation of 
entropy around the equilibrium state and solving the zero-eigenvalue problem, 
the criterion for the onset of gravothermal instability is obtained. The main 
results of our previous analysis are summarized as follows: 

(i) Local entropy extremum ceases to exist in cases with polytrope index 
n > 5 for sufficiently larger radius of the wall, r e > \ crit GM 2 /(—E), and for 
highly density contrast, p c / ' p e > D crit , where M and E denote the total mass 
and energy of the system, p c and p e mean the density at center and edge, 
respectively. 

(ii) The critical values A cr i t and D crit depend on the polytrope index, both of 
which respectively approach 0.335 and 709 in the limit of n — > oo, consistent 
with the well-known result adopting the Boltzmann-Gibbs entropy. 

(iii) The stability/instability criterion obtained from the second variation of 
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Tsallis entropy exactly matches with the result from standard turning-point 
analysis. 

While the successful results suggest that non-extensive generalization of ther- 
modynamics will offer various astrophysical applications involving long-range 
nature of self-gravitating systems, there still remain some important issues 
concerning the physical interpretation of thermodynamic instability. 

Heuristically, the gravothermal instability is explained by the presence of neg- 
ative specific heat as follows. In a fully relaxed gravitating system with suf- 
ficiently larger radius, negative specific heat arises at the inner part of the 
system and we have CV.mncr < 0, while the specific heat at the outer part 
remains positive, Cv.outcr > 0, since one can safely neglect the effect of self- 
gravity. In this situation, if a tiny heat flow is momentarily supplied from 
inner to outer part, both the inner and the outer parts get hotter after the hy- 
drostatic readjustment. Now imagine the case, Cv,outcr > |Cv,inner|- The outer 
part has so much thermal inertia that it cannot heat up as fast as the inner 
part, and thereby the temperature difference between inner and outer parts 
increases. As a consequence, the heat flow never stops, leading to a catastrophe 
temperature growth. 

While the above thought experiment is naive in a sense that we artificially di- 
vide the system into the inner and the outer part, the argument turns out 
to capture an essence of the thermodynamic instability in cases with the 
Boltzmann-Gibbs entropy. Evaluating the specific heat explicitly, Lynden-Bell 
and Wood [7] showed that the specific heat of the total system should be greater 
than zero at the onset of instability, although the central part of this system 
still has the negative specific heat. Therefore, one can naively expect that the 
self-gravitating system generally exhibits the thermodynamic instability asso- 
ciated with the negative specific heat and this could even hold in the system 
characterized by the non-extensive entropy. 

To address this issue, however, we should remember the following two remarks 
that have been never clarified. First note that there exists a subtle point 
concerning the concept of temperature in the non-extensive thermodynamics. 
Framework of the non-extensive formalism is formally constructed keeping 
the standard result of thermodynamic relations [16-18], however, the physical 
temperature, T p h ys , might not be simply related to the usual one, i.e, the 
inverse of Lagrange multiplier, as has been criticized recently [14, 15]. This point 
is in particular important in evaluating the specific heat. 

Second, as has been mentioned by the pioneer work of Lynden-Bell & Wood[7], 
self-gravitating system shows various types of thermodynamic instability. While 
our early study deals with the stellar system confined within an adiabatic wall, 
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one may replace the adiabatic wall with the thermally conducting wall sur- 
rounded by a heat bath. In this situation, assuming the Boltzmann-Gibbs 
entropy, Lynden-Bell & Wood showed that no equilibrium state exists for suf- 
ficiently low temperature and high-density contrast. Note that even in this 
case, the presence of negative specific heat plays an essential role for the ap- 
pearance of instability. 

Keeping the above remarks in mind, in this paper, we focus on the thermody- 
namic property of self-gravitating systems characterized by Tsallis' generalized 
entropy. For this purpose, we first investigate the thermodynamic temperature 
of the self-gravitating system from the Clausius relation. To clarify the physi- 
cal interpretation of thermodynamic instability, the specific heat is computed 
and a role of negative specific is discussed in detail. Then we turn to focus on 
the thermodynamic instability in a system surrounded by the heat bath. The 
stability/instability criterion is derived from the second variation of free energy 
and a geometrical construction of marginal stability condition is discussed. 

While the problem considered here includes some general issues that are com- 
monly faced with the application of the non-extensive thermostatistics, we 
will tackle this problem based on the old Tsallis formalism using the stan- 
dard statistical average, which is currently un-common (e.g., [14,15] [18]). The 
reason why we do not adopt the standard Tsallis formalism using the nor- 
malized g-expectation values is twofold. As has been mentioned in previous 
paper, a naive application of the new formalism apparently shows a problem- 
atic difficulty in our case of the maximum entropy principle (Sec. 2), while no 
such difficulties arise when we apply the earlier formalism (see ref.[2] in de- 
tails). Another reason is that while the new formalism has been deliberately 
constructed so as to eliminate the undesirable divergences in some physical 
systems [18] especially with fractal nature, no serious divergences have ap- 
peared in our case. Precisely speaking, the physical quantities, e.g. mass and 
energy, may have divergence for some equilibrium configuration of the self- 
gravitating system. In order to remedy this divergence, we confine the system 
within a spherical wall, which is a standard prescription in studies for the 
self-gravitating system [6,7]. Since even the old Tsallis formalism preserves a 
consistent framework that recovers the usual thermodynamic structure, from 
a more general view of the non-extensive thermostatistics, we expect that the 
present analysis still provides a valuable insight to the thermodynamic stabil- 
ity of stellar self-gravitating system. Of course, the analysis using normalized 
g-value must play an important role in the Tsallis' non-extensive framework 
and we plan to extend our analysis to the one with the new formulation near 
future. In this sense, present work can be regarded as a preliminary analysis 
toward the next step. This point will be discussed in the last part of this paper, 
together with some implications. 

This paper is organized as follows, in section 2, we recast the problem that finds 
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the most probable state of equilibrium stellar distribution adopting the Tsallis 
entropy. The main part of this paper is section 3, in which the thermodynamic 
properties of stellar polytrope are investigated in detail. After identification 
of the thermodynamic temperature, the explicit expression for specific heat 
is presented and the marginally stability condition for the thermodynamic 
instability is investigated in both the adiabatic and the isothermal cases. In 
section 4, thermodynamic instability in a system surrounded by a thermal 
bath is re-considered by means of the free energy and the marginal stability 
condition is re-derived from the second variation of free energy. Furthermore, 
following the preceding results, the origin of the non-extensive nature in stellar 
polytropic system is discussed in section 5. Finally, section 6 is devoted to the 
summary and discussions. 



2 Stellar polytrope as an extremum state of Tsallis entropy 

In this section, we recast the problem finding the most probable state of equi- 
librium stellar system, based on the maximum entropy principle. In our pre- 
vious study, the entropy for the phase-space distribution function has been 
introduced without recourse to the correct dimensions. Although this does 
not alter the stability/instability criterion for the stellar equilibrium state, 
for the sake of the completeness and the later analysis, we repeat the same 
calculation as shown in ref. [2], taking fully account of the correct dimensions. 

Suppose a system containing N particles which are confined within a hard 
sphere of radius r e . For simplicity, each particle is assumed to have the same 
mass m and interacts via Newton gravity. The problem considered here is to 
find an equilibrium state in an adiabatic treatment. That is, we investigate 
the equilibrium particle distribution in which the particles elastically bounce 
from the wall, keeping the energy E and the total mass M(— Nm ) constant. 

For present purpose, it is better to employ the mean-field treatment that 
the correlation between particles is smeared out and the system can be fully 
characterized by the one-particle distribution function, f(x, v), defined in six- 
dimensional phase-space (x,v) [2,3] [6-9]. Let us denote the phase-space ele- 
ment as h 3 (= IqVq) with unit length l and unit velocity v . Since the distri- 
bution function f(x, v) counts the number of particles in a unit cell of phase- 
space, using the standard definition of the statistical average, the energy and 
the total mass are respectively expressed as follows: 



E = K + U = m 



I 
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M 



m N = m J f{x,v) d % r , 



with the quantity $ being the gravitational potential: 



$(jc) 



-Gm 



•J 



\x — x'\ 



d b r'. 



(2) 



(3) 



In the above expressions, the dimensionless integral measure d 6 T is introduced: 

(4) 



6 d 3 xd 3 v 
d r = — — — ; n = l vo. 



h 3 



Owing to the maximum entropy principle, we explore the most probable state 
maximizing the entropy. The entropy quoted here is a quantity defined in 
the phase-space and it counts the number of possible particle state. We are 
specifically concerned with the equilibrium state for the Tsallis entropy [1]: 



N 



-J 



N 



d b T. 



(5) 



Maximizing the entropy S q under the constraints reduces to the following 
mathematical problem using Lagrange multipliers a and f3: 



5S q -a5M - (35E = 0, 
which leads to [2,12,13]: 



f(x,v) 



1 



1/(9-1) 



(6) 



(7) 



where the constants A and $ are respectively given by 

^ 1/(9-1) 



A = N \ (?—^)m {3 



$ 



1 — (q — l)m a 
(q-l)m o p 



(8) 



The one-particle distribution function (7) is often called stellar polytrope, 
which satisfies the polytropic equation of state [3] [12]. The density profile 
p(r) and the isotropic pressure P(r) at the radius r = \x\ are respectively 
given by 



p(r) = m J f(x,v) 



d 3 v 
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= V27T B (?, ^ {$o - $( r ) } V(«-i) + 3/ 2 , (9 ) 

and 



P(» = ///,, / - fix. v) ^ 



l ° I \ v2 f( x ' v ) 

i ° 5 W 



+ o P(r) {$ -$(r)}, (10) 



\q- 1 2 

with -B(a, 6) being the /3 function. Thus, these two equations lead to the rela- 
tion 

P(r)=K n p 1+1 / n (r), (11) 



with the polytrope index given by 

71=—*— + -. (12) 
q-1 2 V ; 



In equation (11), the dimensional constant K n is introduced: 
1 l\m A V ^ 



n + 



Note that the above quantity is equivalent to the variable (n — 3/2)T/ (n + 1) 
defined in ref.[2]. 

Once provided the distribution function, the equilibrium configuration can be 
completely specified by solving the Poisson equation. Hereafter, we specifically 
restrict our attention to the spherically symmetric configuration for q > l(or 
n > 3/2), in which the dynamically stable state is safely attainable and the 
thermodynamic arguments turn out to capture the physical relevance [3]. 



From the gravitational potential (3), it reads 



Combining (14) with (9), we obtain the ordinary differential equation for $. 
Equivalently, a set of equations which represent the hydrostatic equilibrium 
are derived using (9), (10) and (14): 
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dP(r) Gm(r) 
dm(r) 



dr 



47rp(r) r 2 . (16) 



The quantity m(r) denotes the mass evaluated at the radius r inside the wall. 
We then introduce the dimensionless quantities: 



which yields the following ordinary differential equation: 

9" + jff + 9 n = 0, (18) 



where prime denotes the derivative with respect to £. The quantities p c and P c 
in (17) are the density and the pressure at r = 0, respectively. To obtain the 
physically relevant solution of (18), we put the following boundary condition: 

0(0) = 1, 0'(O)=O. (19) 



A family of solutions satisfying (19) is referred to as the Emden solution, which 
is well-known in the subject of stellar structure (e.g., see Chap. IV of ref.[19]). 

Figure 1 shows the numerical solution of equation (18) for various polytrope 
indices, where the density profile, p{r)/p c is plotted as a function of dimen- 
sionless radius, £. Clearly, profiles with index n < 5 rapidly fall off and they 
abruptly terminate at finite radius(/e/t-pcme/), while the n > 5 cases infinitely 
continue to extend over the outer radius(right-panel) . As already mentioned 
in previous paper, characteristic feature seen in figure 1 plays an essential role 
for the thermodynamic instability associated with negative specific heat. 

For later analysis, it is convenient to introduce the following set of variables, 
referred to as homology invariants [19,20]: 



_d In m(r) _ 4nr 3 p(r) _ £8 n 
d\nr m{r) 9' ' 

_ dlnP(r) p(r) Gm{r) = _ (n + 1) ^ (21) 



dlnr P(r) r 9 
which reduce the degree of equation (18) from two to one. The derivative of 
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these variables with respect to £ becomes 

du / n \ u dv ( „ 1 \ f 
_=(3_ u u] _ = (_i + u+ 22 

rf£V n+1 7 f ' rf£V n+1 7 £ v; 

Equations (18) can thus be re- written with 

u dv_ = (n + l)(u- l)+v ^ 
v du (n + 1)(3 — u) — nv 

The corresponding boundary condition to (19) becomes (u,v) = (3,0). Using 
these variables, the basic thermodynamic quantities such as the energy and 
the entropy are evaluated and the results are summed up in Appendix A, 
which are subsequently used in section 3. 



3 Thermodynamic properties of stellar polytrope 

In this section, we address our main issue, i.e, the physical interpretation of 
gravothermal instability in stellar polytropes, based on the framework of non- 
extensive thermodynamics. In section 3.1, we first discuss the thermodynamic 
temperature of stellar polytrope calculating both the heat and the entropy 
changes in a quasi-static treatment. Then we evaluate the specific heat in 
section 3.2. The connection between the absence of extremum entropy state 
and the presence of negative specific heat is discussed in detail. Further, we 
argue that there appears another type of thermodynamic instability, which is 
subsequently analyzed by means of the free energy. 

3. 1 Thermodynamic temperature from the Clausius relation 

As has been mentioned in section 1, the concept of temperature is non-trivial 
in non-extensive thermostatistics. This is because the standard framework of 
thermodynamics crucially depends on the assumption of extensivity of entropy. 
According to the recent claim, the definition of physical temperature T p h ys 
should be altered depending on the choice of energy constraint and is related 
to the inverse of the Lagrange multiplier, with some correction factors 
[14,15]. Note, however, that this discussion heavily relies on the extensivity 
of the energy as well as the thermodynamic zeroth law. In our present case, 
the maximum entropy principle was applied subject to the constraints E and 
M, adopting the standard definition of mean values (see eqs.(l)(2)). As a 
consequence, the resultant energy E becomes non-extensive and we cannot 
apply the above definition. 
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To address the physical temperature in the present case, we therefore consider 
the relation between the heat transfer and entropy change and seek the most 
plausible candidate for thermodynamic temperature. That is, we analyze the 
variation of equilibrium configuration under fixing the total mass. Specifically, 
we deal with the quasi-static variation along an equilibrium sequence. 

Let us first write down the heat change. The thermodynamic first law states 
that 

d'Q = dE + P e dV, (24) 



where the operation d! stands for incomplete differentiation. The subscript e 
denotes a quantity evaluated at the edge. In the spherically symmetric con- 
figuration, the second term in right-hand side of (24) becomes ^r 2 e P e dr e . As 
for the first term, the energy of the stellar polytropic system within the radius 
r e , is computed in Appendix A.l. Introducing the dimensionless parameter A, 
it is expressed in terms of the homology invariants as follows: 



A: 



GM 2 



n 



!{l-( B+ l)I} + (n -2)£ 



(25) 



where the quantity with subscript e represents the one evaluated at the bound- 
ary r = r e . Using (25), the heat change d'Q is rewritten as follows: 



/ CM 2 \ 
d'Q = d[ -A + 4vrr e 2 P e dr e , 



GM 2 |A + «e\ dr e dX df e j 

where the relation Airr^P e /(GM 2 ) = u e /v e is used in the last line (see defini- 
tions (20)(21)). In the above expression, derivative of A with respect to £ e can 
be computed with a help of relation (22) (see eq.(33) of ref.[2]): 

d\_ n-2 g(u e , v e ) 
a£ e n — 5 2v e 



where 

gM = + 2», - { 8 + 3 (=±1)}„- ^„ + 3 (2±|) . (28) 



Next focus on the change of the entropy. From (71) in Appendix A. 2, the 
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entropy of the extremum state is given by 



S„— [ n - - 



1 PGM 2 



n 



[ v e v e J 



(29) 



Hence, the variation of entropy dS q under fixing the total mass can be decom- 
posed into the variation of homology invariants (u e , v e ), radius r e and Lagrange 
multiplier j3 as follows: 



dS n 



n-3/2 PGM 2 



n — 5 



(f-ti) ; 2 «i-( n+ i)i + i 

V P r e J { v e v e 

(-Ue / dUe_ _ dVe\ _ U + 1 dv e 
\ V e \U e V e j V e V e 



(30) 



Among these variations, variation of homology invariants is simply rewrit- 
ten with d£ e , through the relation (22). On the other hand, from the mass 
conservation, the variation of Lagrange multiplier, dp is related to both the 
variations of homology invariants and dr e as follows. Using the condition of 
hydrostatic equilibrium at the edge r e , one can obtain the following relation 
(see derivation in Appendix A. 3): 



_ \ {GMY{ mo PY-*l 2 \ 1,(n ~ l) ^ 
V= { J =a n (u e v e ) /( >, 



(31) 



where the constant a n is given by 

f (n-l/2)«- 3 / 2 



l/(n-l) 



a. 



16v^tt 2 (n + l) n 5(3/2, n - 1/2) 



(32) 



which asymptotically approaches unity, in the limit n 
total mass M constant, variation of (31) yields 



n — 3/2 dp n — 3 dr e 
n — 1 P n — 1 r e n — 1 \ u 



1 / du e dv e 
I h n — 

V P 



+oo. Keeping the 



(33) 



We then rewrite it with 
1 



dp 
J 



dr e 

r P . 



n 



3/2 



3 dr e du e 

2 7V Ue 



+ n 



dVe 



(34) 



Substituting the relation (34) into equation (30), the dependence of dp/P can 
be eliminated. Thus, using the relation (22), the final form of the entropy 
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change is expressed in terms of the variations d£ e and dr e . After some manip- 
ulation, we obtain 



dS„ 



(3GM 2 



3/2 



n — 5 



u e n + 1 \ dr e n — 2 1 
2— + 11 — 

v e v e J r e n — 5 2v e 



x Uut 



2u e v e 



8 + 3 



n+1 



n 



v e + 3 



n 



n+ 1 



n 



• (35) 



Now, from the knowledge of the expressions A and £ e (d\/d£ e ), one can easily 
show that the above equation is just identical to 

dSq = ^Ml{( X + ^)±l-^f\. (36) 
r e [V v e J r e d£ e £ e J 



Therefore, comparison between (36) and (26) immediately leads to the follow- 
ing relation: 

dS q = [3d'Q = !3 (dE + P e dV), (37) 



which exactly coincides with the standard result of Clausius relation in a 
quasi-static process. 

The relation (37) strongly suggests that the thermodynamic temperature T p h ys 
is identified with the inverse of Lagrange multiplier, T p h ys = 1//3- At first 
glance, the result seems somewhat trivial, since one can easily expect this rela- 
tion from the standard thermodynamic relation, 8S q /dE = (3, which generally 
holds even in the the non-extensive Tsallis formalism [16,17]. As advocated by 
many author, however, the relation OS q /dE = f3 does not simply imply the 
thermodynamic temperature T p h ys = 1//3 and it might even contradict with 
the thermodynamic temperature defined through the thermodynamic zeroth 
low [15]. 

On the other hand, in our case of the self-gravitating system, the thermody- 
namic temperature T phys = 1/(3 is mathematically verified by the integrable 
condition of the thermodynamic entropy through the Clausius relation. Fur- 
ther, it is remarkably found that the relation T p h ys = 1/(3 holds even in the 
absence of gravity (the limit G — * 0) and can be proven through an alterna- 
tive route. In Appendix B, as a pedagogical example, we demonstrate that the 
relation T phys = 1/(3 is indeed obtained in the classical gas model using the 
Carnot cycle. 
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3.2 Negative specific heat and thermodynamic instability 



Once obtained the thermodynamic temperature, T phys = 1//3, we are in a 
position to investigate the thermodynamic instability from the straightforward 
calculation of the specific heat. Let us first discuss the qualitative behavior of 
the specific heat. By definition, the specific heat at constant volume is given 
by 




Recall that the dimensionless parameters A and r] are respectively proportional 
to —E and /jC"- 1 )/^ 3 / 2 ) ( see eqs.(25)(31)). This implies that for a system of 
constant mass inside a fixed wall, the qualitative behavior of (38) can be 
deduced from the relation between 77 and A. 

Figure 2 depicts the trajectories of the Emden solutions in the (77, A)-plane 
with various polytrope indices. Each point along the trajectory represents 
an Emden solution for different value of the radius r e . From the boundary 
condition, all the trajectories start from (77, A) = (0, — 00), corresponding to 
the origin r e = 0. As gradually increasing the radius, the trajectories first 
move to upper-right direction monotonically, as marked by the arrow. At this 
stage, the kinematic energy dominates the potential energy and the system 
lies in a kinematically thermal state (A < 0), indicating the positive specific 
heat. For larger radius, while the curves with index n < 3 abruptly terminate, 
the trajectories with n > 3 suddenly change their direction from upper-right 
to upper-left. Moreover, in the case of n > 5, the trajectory progressively 
changes its direction and it finally spirals around a fixed point. 

From these observations, one can roughly infer the existence of the two types 
of the thermodynamic instability as follows. At first inflection point for n > 3, 
the specific heat diverges and the signature of C v becomes indefinite. Beyond 
this point, the specific heat changes from positive to negative. This means 
that the potential energy conversely dominates the kinetic energy, indicating 
the system being gravothermal. In this case, equilibrium state ceases to exist 
for a system in contact with a heat bath, but does still exist for a system 
surrounded by an adiabatic wall. However, for the polytrope index n > 5, the 
specific heat of the system turns to increase beyond this inflection point and 
it next reaches at the point dX/dr] = 0, i.e, C v = 0. This means that while 
the inner part of the system still keeps the specific heat negative, the fraction 
of the outer normal part grows up as increasing r e and it eventually balances 
with inner gravothermal part. Thus, beyond this critical point, no thermal 
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balance is attainable and the system becomes grave-thermally unstable. This 
is true even in the system surrounded by an adiabatic wall. 



Now, let us write down the explicit expression for the specific heat C v . In 
equation (38), the variation of f3 and E with £ e can be respectively rewritten 
with 

and 

(d$\ n-l g dq 
{d£)~ n-3/2^ 



(40) 



Here, the variable d\/d^ e has been already given in (27). As for the derivative 
of rj with respect to £ e , we obtain 

e. ? = U - n -^) * (4i) 



£ e V n-l 
Then the quantity C v becomes 

(n-3/2)(n-2) [3GM 2 g(u e ,v e ) 



C v = 



(n-l)(n-5) r e , / _ n^3 



2v e M e 

n 



with the function g(u e ,v e ) given by (28). Notice that the above expression is 
still redundant, since there remains the explicit dependence of the variable j3. 
Eliminating the variable (5 by using the relation (31), one finally obtains 

CV _ . (3/2)/( "~ 3/2) 9M) , v3/2) V(n-3/ 2 ) 

iV GMrJ n-3\ \ Ue v - ) ' ^ 



2 |u e - 



n 



where we introduced the new dimensionless constant a n : 

a„ S ( ';-^' ( "-f aV<-V* (43) 
(ra — — 5) 

Note that in the limit n — > +oo, equation (42) consistently recovers the well- 
known result of isothermal sphere (e.g, eq.(39) of ref.[22]): 

n— >+oo 4w e + 2u e v e — ll"U e + 3 / a a\ 

~W ~" 2(u e - 1) • 1 j 
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Comparing (42) with the isothermal limit, the resultant expression contains a 
residual dimensional parameter h, as well as the quantities M and r e . While 
the residual dependence can be regarded as a natural consequence of the non- 
extensive generalization of the entropy, it would be helpful to understand the 
origin of this scaling in more simplified manner. This will be discussed in 
section 5. 

Apart from the residual factor, the expression of specific heat (42) clearly 
reveals the two types of thermodynamic instability seen in Figure 2. The 
inflection point with the infinite specific heat, C v — > ±oo leads to the condition 

u e - ^—^ = 0, (45) 
n — 1 



which immediately yields the conclusion that this is only possible for the 
polytrope index n > 3, consistent with Figure 2. On the other hand, critical 
point with the vanishing specific heat, C v = corresponds to the following 
condition: 

g(u e , v e ) = 0. (46) 



This is exactly the same condition as obtained from the second variation of 
entropy (see eq.(33) or (53) in ref.[2]). According to the previous analysis, the 
condition (46) represents the marginal stability at which the extremum state 
of the entropy S q is neither maximum nor minimum. This situation turns out 
to appear when the polytrope index n > 5. 

Therefore, we reach a fully satisfactory conclusion that the thermodynamic in- 
stability found from the second variation of entropy is intimately related to the 
presence of negative specific heat and the stability/instability criterion can be 
exactly recovered from the critical point of the thermal balance, C v = 0, which 
is also consistent with the analysis in the Boltzmann-Gibbs limit, n — > oo [7]. 
The successful result can be regarded as an outcome of the correct definition 
of T phys . As for the transition point with C v — > ±oo, it clearly indicates the 
thermodynamic instability of a system in contact with a thermal bath. In next 
section, by means of the free energy, we confirm that the condition (45) indeed 
represents the marginal stability of the system surrounded by a thermal wall 
and beyond this point the system will be unstable. 

In Figure 3, by varying the radius r e , the normalized specific heat per particle 
C*/N is plotted as a function of density contrast, p c /p e around the critical 
polytrope indices n = 3(upper-panels) and n = ^(middle-panels) . Here, the 
normalized specific heat C* is defined by the specific heat C Y divided by 
the redundant factor (/i 2 /GMr e )^ 2 '^" _3 / 2 '. Obviously, the transition point 
C v — > ±oo appears when n > 3(crosses), while the existence of critical point 
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C v = is allowed for higher density contrast of n > 5 cases (arrows). The 
critical values -D C rit = (Pc/Pe)crit indicated by arrows exactly coincide with 
those obtained from the previous analysis (see Table 1 of ref.[2]). Lower-panels 
of Figure 3 show the specific heat with large polytrope indices n — 10 and 
30, together with the Boltzmann-Gibbs limit (n — > +oo, labeled by iso). As 
increasing the polytrope index n, the critical/transition points tend to shift 
to the lower density contrast, while the successive divergent and zero-crossing 
points appear at the higher density contrast, corresponding to the behavior 
seen in Figure 2. 



Thermodynamic instability from the second variation of free en- 
ergy 



Previous section reveals that there exists another type of thermodynamic in- 
stability in which the marginal stability is deduced from the condition (45). 
In this section, to check the consistency of the non-extensive thermostatistics, 
we reconsider this issue by means of the Helmholtz free energy: 

F q = E — T phys S q . (47) 



Adopting the relation T p h ys = 1//3, we re-derive the marginal stability condi- 
tion (45) from the second variation of F q . 

Consider a system surrounded by the thermally conducting wall in contact 
with a heat bath. Usually, the stable equilibrium state should keep the free 
energy F q minimum. Thus the presence of thermodynamic instability implies 
the absence of minimum free energy, which can be deduced from the signature 
of the second variation 5 2 F q around the extremum state of free energy. Since 
the non-extensive formalism still verifies the Legendre transform structure 
leading to the standard result of thermodynamic relation[16,17], the extremum 
state of the free energy exactly coincides with that of the entropy. One thus 
skips to find the extremum state of F q and proceeds to evaluate the second 
order variation. 

In contrast to the adiabatic treatment, we here deal with the density pertur- 
bation p — > p + 5p, surrounded by a thermal wall. To be specific, we evaluate 
the second variation under keeping the radius r e , the total mass M and the 
temperature T p h ys constant. Then the variation of energy up to the second 
order leads to 



5E= 5 



J {-P(x) + -p(x)$(x)}d 3 x 
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7 {l8P+\{6p* + p6*) + ±6 P 6*}<? 



X. 



(48) 



Similarly, using the expression (70) in Appendix A. 2, the variation of Tsallis 
entropy becomes 



6S q = 5 



n — 



— — \ n 



N-/3 J P{x)d 3 x 
(3 J 5P(x)d 3 x. 



(49) 



The above expressions include the variation of pressure 5P, which can be 
expanded with a help of the polytropic equation of state (11): 



5P 



1 + 



1\ P 



nj p 



* P+ i(i + 



nj n p A 



(50) 



Combining the above result with equations (48) and (49) and collecting the 
second order terms only, the second variation of free energy becomes 



8 2 F q = 5 2 E - T phys 5 2 S n = ^ 



1 

q ~ 2 



n+1 P 



n p* 



- (Spy + <5p5$ d s x, (51) 



where the relation T p h ys = is used in the last line. Now, restricting our 
attention to the spherical symmetric perturbation, we introduce the following 
perturbed quantity (see refs. [2] [8]): 



Sp(r) 



1 dQ{r) 



4-7rr 2 dr 



(52) 



Then the mass conservation SM = implies the boundary condition Q(0) = 
Q(r e ) = 0. Substituting (52) into (51) and repeating the integration by part, 
one finally reaches the following quadratic form: 



1 r ' e 

5 2 F q = -- fdr Q(r) 



n+1 d 



P\ d 



G 



n 



dr { 4-7T r 2 p \p J dr 



Q(r). (53) 



Thus, the problem just reduces to the eigenvalue problem and the stability 
of the system can be deduced from the signature of the eigenvalue. More 
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specifically, the onset of instability corresponds to the marginally stability 
condition, 5 2 F q = 0, and it is sufficient to analyze the zero-eigenvalue equation: 



L Q(r) = 



d r i 

dr \ Anr 2 p 




Q{r) = 0, 



(54) 



with the boundary condition, Q(0) = Q(r e ) = 0. Equation (54) has quite sim- 
ilar form to the zero-eigenvalue equation found in the adiabatic treatment (see 
eq.(46) of ref.[2]). Except for the non-local term, one can utilize the previous 
knowledge to solve the equation (54): 

f 3 n-3Gm(r) f n — lGm(r) 

L 47rr 3 p) = — — L m r = — — 55 

n + 1 r l n + 1 r l 

These two equation leads to the ansatz of the solution: 

f 7i — 3 1 
Q(r)=c< 47rr 3 p(r) m(r) > . (56) 



Here, the variable c is an arbitrary constant. The above equation (56) au- 
tomatically satisfies the boundary condition Q(0) = 0, while the remaining 
condition Q(r e ) = puts the following constraint: 

Q(r e ) = c Unrl p e - — ^ m]=c(u.- — | ) M = 0. (57) 
V n — 1 / V n — 1 / 



Again, we arrive at the satisfactory result that the solution of zero-eigenvalue 
equation exactly recovers the condition (45). 

Now, remaining task is to show that the second variation 5 2 F q becomes nega- 
tive beyond the transition point of C v — > ±oo. One can rewrite the expression 
(53) with 




o 

with the constant H given by 




o 
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That is, the condition H > 1 implies stable local minimum state of free en- 
ergy, while the inequality H < 1 represents unstable local maximum state. 
Integrating by part, equation (58) can be regarded as an eigenvalue equation 
with eigenvalue, H: 




H — ^— r (59) 
n + 1 r 2 



Obviously, equation (56) becomes the solution of above equation with the 
minimum eigenvalue, H m \ n = 1, if the condition (57) is fulfilled. In this case, 
solution (56) can be regarded as the ground state of the eigensystem (59), 
since the function (56) does not possess any nodes between [0,r e ]. Therefore, 
for a suitably smaller radius r e or a smaller density contrast p e j p c below the 
transition point, the eigenvalue H should be larger than unity. Conversely, 
from continuity, the condition H < 1 must be satisfied beyond the critical 
radius. 

Finally, using the (u, v )-variables, the geometrical meaning of onset of ther- 
modynamic instability is briefly discussed in similar manner to the adiabatic 
case. In Figure 4, the thick solid lines show the Emden trajectories with various 
polytrope indices in (u, v )-plane. The thin-solid lines in Figure 4 represents the 
straight lines, u — (n — 3)/(n — 1) = 0. Since the equilibrium state only exists 
along the Emden trajectory, the condition (57) is satisfied at the intersection 
of these two solid lines, which is only possible for n > 3. On the other hand, 
as seen in previous section, the equilibrium system surrounded by a thermal 
wall is characterized by the three parameters, r e , M and (3(or T phys ), through 
the relation (31). In other words, the system must lie on the curve: 

n \ (n-l)/n 

■^-i u , (60) 



Or 



with some constant value 77. We have seen in Figure 2 that the constant value 
i] is bounded from above, r\ < r] crit . Thus, the critical curve (60) with r\ = i] crit 
must intersect with both the Emden trajectory and the straight line u — (n — 
3)/(n — 1) = simultaneously. This is clearly shown in Figure 4, where the 
critical curve is plotted as dashed lines. Since the critical curves tangentially 
intersect with Emden solutions, it always satisfies the condition drj/d^ = at 
the contact point, leading to the condition (45) consistently. 

Table 1 summarizes the dimensionless quantities 77 crit and D CTlt = (p c /p e )crit 
evaluated at the contact point. As increasing the polytrope index n, these 
values asymptotically approach the well-known results of Boltzmann-Gibbs 
limit, r] crit — > 2.52 and D CTit — > 32.1. 
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5 Origin of non-extensive nature in stellar polytrope 



As has been mentioned in section 3.2, specific heat of the stellar polytropic 
system explicitly depends on the residual dimensional parameter h, in contrast 
to the isothermal limit (44). In this section, to contact the physical meaning of 
the non-extensivity in stellar polytrope, we discuss the origin of this residual 
dependence. Indeed, the appearance of the residual factor can be recognized 
as the breakdown of both the intensivity of temperature and the extensivity 
of energy and entropy as follows. From equation (18), the asymptotic behavior 
of the Emden solution becomes 



_ ^-2/(n-l) > p „ r -2n/(n-l) > ^ r ^ Qo) 

so that the mass within a sphere of radius r is given by 

M ~pr 3 oc r (»-3)/(»-i) . (61) 

Then the energy of a virialized stellar system is roughly estimated as 



E ~ oc r-t"- 5 )/^- 1 ) oc M (n " 5)/(n - 3) 

and the relation (31) tells 



P OC r -(™-3)/(«-l)/(«-3/2) K M -l/(n-3/2)_ 

These relations clearly show the breakdown of the intensivity of temperature 
and the extensivity of energy, which lead to the scaling of the specific heat per 
mass: 



Cy _ I dE I3E_ w-3(n-2)/(n-3)/(n-3/2) (62) 

N~MdT phys M ■ {bI) 

On the other hand, the dimensionless combination h 2 / (GMr e ) represents the 
ratio of a typical scale of the stellar system, GMr ~ (GM/r)r 2 ~ v 2 r 2 , to 
that of the reference cell, h = v l . This behaves as 



oc _1_ oc M 2(n-2)/(n-3)_ ( g 3) 



GMr e Mr, 
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Thus, these two equations (62) and (63) lead to the scaling relation of (42): 



r / 7,2 \ (3/2)/(n-3/2) 

• (64) 



N \ GMr 



e , 



Notice that the Clausius relation (37) suggests that the entropy per unit mass 
has the same scaling relation: 



S q (3E a 



v 



M M N ' 

Therefore, resultant dependence (64) for the stellar polytrope can be a natural 
outcome of the non-extensivity of the entropy. 

In fact, framework of the thermostatistics generally requires an introduction 
of the scale of the unit cell in order to count the available number of states in 
phase spaces. This is even true in the case of the isothermal stellar system(n — > 
+oo or q — > 1), but, the thermodynamic quantities show somewhat peculiar 
dependence of the scale h. A typical example is the entropy: 



BG m \\ e e 2J \ 4n J 2 ^TiGMrej] 1 

where u e and v e are the homology invariants for the isothermal system. The 
above equation shows that in the Boltzmann-Gibbs limit, /Vdependence of the 
entropy can be recognized as a matter of choice of an additive constant, so 
that its derivatives, e.g., specific heat, is free from the residual dependence. 

It should be emphasized that the stellar equilibrium system recovers the ex- 
tensivity in the limit n — > oo and it behaves as 



E ~ M ~ r, Cy ~ M. (65) 

Also, the temperature becomes intensive in this limit. Thus, we readily under- 
stand that the scaling behavior shown in (42) or (64) has nothing to do with 
the long-range nature of the gravity. Even in the free polytropic gas model in 
Appendix B, the residual dependence emerges as 



~W ~ 1 (p/ p )v 2 / 3 j 

It follows that the explicit dependence of the specific heat on the reference cell 
scale h just originates from the the non-extensive nature of Tsallis entropy. 
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6 Summary &: Discussions 



In this paper, thermodynamic properties of the stellar self-gravitating sys- 
tem arising from Tsallis' non-extensive entropy have been studied in detail. 
In particular, physical interpretation of the thermodynamic instability previ- 
ously found from the second variation of entropy is discussed in detail within 
a framework of the non-extensive thermostatistics. After briefly reviewing the 
equilibrium state of Tsallis entropy, we first address the issues on thermo- 
dynamic temperature in the case of equilibrium stellar polytrope. Analyzing 
the heat transfer and the entropy change in a quasi-static process, standard 
form of the Clausius relation is derived, irrespective of the non-extensivity of 
entropy. According to this result, we explicitly calculate the specific heat and 
confirm the presence of negative specific heat. The onset of instability found in 
previous work just corresponds to the zero-crossing point, C v = 0, supporting 
the fact that the heuristic explanation of gravothermal catastrophe holds even 
in the non-extensive thermostatistics. 

Further, the analysis of specific heat shows divergent behavior at n > 3, 
suggesting another type of thermodynamic instability, which occurs when the 
system is surrounded by a thermal wall. We then turn to the stability analysis 
by means of the Helmholtz free energy. Similar to the previous early work, the 
stability/instability criterion just reduces to the solution of the zero-eigenvalue 
problem and solving the eigenvalue equation, we recover the marginal stability 
condition derived from the divergence of specific heat (45). 

In addition to the thermostatistic treatment, we have also discussed the ori- 
gin of non-extensivity in stellar polytrope. The residual dependence of the 
reference scale h appeared in the specific heat (42) naturally arises from the 
non-extensivity of the entropy and the resultant scaling dependence can be 
simply deduced from the asymptotic behavior of the Emden solutions. 

The stability analysis using the free energy in section 4 is consistent with re- 
cent claim by Chavanis [23], who has investigated the dynamical instability of 
polytropic gas sphere. According to his early paper [22], the thermodynamic 
stability of stellar system is intimately related to the dynamical stability of 
gaseous system, which has been clearly shown in the case of the isothermal 
distribution. Thus, the correspondence between Chavanis' recent result [23] 
and a part of our present analysis can be regarded as a generalization of his 
early work to the polytropic system. Note, however, that starting from the 
Tsallis entropy, we extensively discuss the thermodynamic temperature and 
the specific heat of stellar polytrope. Therefore, at least, from the thermody- 
namic point of view, our present analysis provides a valuable insight to the 
stellar equilibrium systems. 
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A particular interest in the thermodynamic relation is the Clausius relation 
(37) that has been still preserved in the non-extensive stellar system. This 
is indeed consistent with the standard thermodynamic relation dS/dE = f3, 
if one keeps the volume constant. Note also that the relation dS/dE = (3 is 
readily obtained from the standard Legendre transform structure. While we 
only dealt with a specific case with the non-extensive entropy (5), it is well 
known that the standard Legendre transform structure does generally hold 
independently of the functional form of the entropy [17]. Hence, our result in 
turn suggests that the Clausius relation is also valid for any stellar system 
maximizing the entropic functional more general than Tsallis'. 

At present, the results shown in this paper seems fully consistent with the 
general framework of the thermostatics. Apart from the thermodynamic in- 
stability, the stellar polytropic system can be a plausible thermodynamic equi- 
librium state, as well as the isothermal stellar distribution. In the isothermal 
case, existence of the thermodynamic limit has been discussed by de Vega and 
Sanchez [21]: 

M 

M, V -> oo, = fixed, 

where V ~ r 3 is a volume of the system. Recalling the discussion in section 5, 
the above condition merely reflects the extensivity of the isothermal system 
(65). Thus, similar argument can hold for the non-extensive system. According 
to the scaling relation (61), the existence of the thermodynamic limit in stellar 
polytrope yields the condition: 

M 

M,V-^oo, v(n _ 3)/{3n _ 3) = fixed. 

Note, however, that this discussion relies on the non-uniqueness of the Boltzmann- 
Gibbs theory, which can be proven only mathematically [24]. Indeed, frame- 
work of the thermostatistics cannot answer the question whether the stellar 
polytropic distribution is really achieved as a thermodynamic equilibrium. 
To address this issue, we must study the detailed process of the long-term 
stellar dynamical evolution. In the light of this, the analysis using Fokker- 
Planck model or direct N-body simulation can provide an invaluable insight 
to the non-extensive nature of stellar gravitating systems. This issue is now 
in progress and will be presented elsewhere. 

Another remaining issue is the re-examination of the present analysis from a 
view of the 'standard' Tsallis formalism using the normalized g-expectation 
values. Apart from some technical issues on the treatment of the maximum 
entropy principle, one might naively expect that the consistency between the 
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statistical and the thermodynamic analysis should be preserved even in the 
new formalism. However, a rather subtle point would be the identification 
of the thermodynamic temperature. As several author stated, the standard 
Clausius relation should be modified in the new Tsallis formalism and the 
resultant form of the expression apparently seems to contradict with the ther- 
modynamic temperature defined through the thermodynamic zeroth law [15]. 
This point will be in particular important in discussing the thermodynamic 
instability and should be clarified along the line of our present treatment. 

This work is supported in part by the Grant-in- Aid for Scientific Research of 
Japan Society for the Promotion of Science. 
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Appendix A: Thermodynamic variables in a stellar polytropic sys- 
tem 



In this appendix, using the equilibrium state of stellar polytrope described 
in section 2, we explicitly evaluate the thermodynamic variables, which have 
been used in section 3 and 4. 



A.l Energy 



Recall that the equilibrium system confined in a spherical container satisfies 
the following virial theorem (e.g, p. 502 of Ref.[3]): 



2K + U = 4irr 3 e P e . 



The energy (1) is then expressed as 



E = K + U = Anrl P e - K = 4nr 3 e P e 



r 2 dr. 



(66) 



To evaluate the above integral in the spherically symmetric case, we use the 
following integral formula: 



jp(r) 



4nr 2 dr = 



If, , MP e GM 2 

\ 8nr 3 e P e - n + 1 - + 

n — 5 I p e r e 



(67) 



which can be derived from the conditions of hydrostatic equilibrium, (15) and 
(16) (see Appendix A of ref. [2]). Thus, the energy of extremum state becomes 



E = 



n — h 



3 f GM 2 . , . MP e 

n+1 e - 

2 { r e p e 



(n - 2) 4:%r 3 e P ( 



In terms of the homology invariants, we obtain 
1 GM 2 



E 



n — 5 r P 



(68) 



(69) 



A. 2 Entropy 



First note the definition of Tsallis entropy (5): 
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*--(»-§)/ 



/ , x (n-l/2)/(n-3/2) 



Substituting the distribution function (7) into the above equation, after some 
manipulation, we obtain 



S q = -[n-~ 



J/3 J P(x) d 3 x - 7vj . (70) 
Thus, the substitution of integral formula (67) immediately leads to 



n — 5 



which can be expressed in terms of the homology invariants: 



S„ 



3 

n 

2 



1 PGM 2 



n — 5 r P 



f 2 ^-(n + l)- + l\ + N 
I v e v e J 



(71) 



Radius-mass-temperature relation 



The mass-radius-temperature relation (31) is derived from the equilibrium 
stellar polytropic configuration. Using (15), we first write down the condition 
of hydrostatic equilibrium at the boundary r e : 



GM 1 (dP\ 
r\ p e \dr) e ' 

The right-hand-side of this equation is rewritten with dimensionless quantities 
in (17): 



GM 



= -(n+l)K n p L e '» (^j 9' e . 



(72) 



We wish to express the above equation only in terms of the variables at the 
edge. To do this, we eliminate the residual dependences, p c and K n from (72). 
The definition (17) leads to 



k 

r P 



AnGp 



1/2 



4ttG 



1/2 



n + 1)P C 



(n + l)K n 



(n-l)/(2n) 
He ; 
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which can be rewritten with 



4nG 



l/(n-l) 



(n+l)K n 



2/(n-l) 



Substituting the above relation into (72), the p c -dependence is first eliminated 
and we obtain 



Qn/(n-l) M 

(n-3)/(n-l) 
I e 



{{n + l)K n y 
An 



l/(n-l) 



^(n+l)/(n-l) 



(73) 



As for independence, the definition (13) together with (8) yields 



<„ + 1) K. = {4V2. *g^> M }-'" (m , r M», (74) 

Hence, substituting the above expression into (73), the relation between mass 
M, radius r e and Lagrange multiplier (3 can be finally obtained. In terms of 
the homology invariants, it follows that 



j(GMnm pr-^\ 1/(n - 1] ^ fii ^ 



\ ' " r^h 3 J 
where the constant a n is given by 



(75) 



«n = 



(n-1/2)"" 3 / 2 1 1/(n_1) 

16V27T 2 (n + l) n 5(3/2, n - 1/2) J 



which asymptotically approaches unity in the limit n — > oo. 



Appendix B: Thermodynamic temperature of classical gas model 
from the Carnot cycle 



In a standard framework of thermodynamics, the temperature is defined by 
means of an efficiency of the Carnot cycle. Here we apply the standard pro- 
cedure to seek the physical temperature T phys for so-called polytropic system 
of which distribution function is given by the extremization of the Tsallis en- 
tropy (see eqs.(5)(6)). For simplicity, we discuss a case of the free classical gas 
without gravity, which corresponds to the G — > limit of the stellar polytropic 
system. 
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From the G — > limit of the formula (68), free polytropic system of the volume 
V with homogeneous pressure P and density p has an (internal) energy: 

E = K = l PV=\—- (™) 
2 2 p y ' 



Here we drop the subscript e for the pressure and density, since both are 
constant within the system in absence of gravity. And equation of state (11) 
becomes 

/M\ 1+1 / n 

P = K n pi+Vn = Kn ^ . (77) 



From equations (8) and (13), the constant K n is related to the Lagrange 
multiplier (5 as 

K n oc /5-(«- 3 / 2 )/™, (78) 



so that this constant can be used as a parameter which characterizes the 
temperature of the system. However, it is not sure whether K n itself has a 
role of the physical temperature, which should be determined through the 
efficiency of the Carnot cycle. 

The internal energy (76) and the equation of state (77) give the thermody- 
namic first law: 



d'Q = dE+ PdV 

= M i + vn J 3 dKn + (^va K jy_\ m 

from which adiabatic changes d'Q = is expressed as 



K n V (2/3 - 1/n) = constant, P V 5/3 = constant'. (80) 

Note that adiabatic lines in a P-V plane become steeper than isothermal ones 
when n > 3/2. 

Now, let us consider the Carnot cycle shown in Figure 5. As usual, quasi- 
static changes B — > C and D — > A are adiabatic. As for the process A — > 
B, the system is in a thermal contact with a heat bath which has a higher 
temperature ■ Similarly, during the change C — > D, the system lies in a 
thermal equilibrium with another heat bath that has a lower temperature K%- 
The system absorbs amount of heat Q H from the higher temperature bath and 
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disposes Q L to the lower one during the isothermal processes B and C- 
D, respectively. They are easily evaluated from (79): 



QH = {n _ |) M l+l/n R H (y-l/n _ y-l/n^ 

Q L — (n— |) M 1+1 /« {y D l/n - V c 1/n ). 



(81) 



On the other hand, a relation between the parameters of the cycle can be 
obtained from the equation of state (77) and the adiabatic changes (80): 



Thus, equations (81) and (82) lead to the following efficiency of the Carnot 
cycle: 

=1 -W< (83) 



where we used the relation (78) in the last line. This clearly shows that the 
inverse of the Lagrange multiplier (5 has a role of the physical temperature. 
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Table 1 

Critical values of the radius-mass-temperature relation, r/ cr it and the density contrast 
between center and edge, -D C rit = (p c /p e )crit i n the case of a system in contact with 
a heat bath for given polytrope index n or q. 



n 


q 


f?crit 


-Dcrit 


3 


5 
3 






4 


7 
5 


0.9421 


153.5 


5 


9 
7 


1.193 


88.15 


6 


1.22 


1.379 


68.38 


7 


1.18 


1.520 


58.86 


8 


1.15 


1.631 


53.28 


9 


1.13 


1.720 


49.62 


10 


1.12 


1.793 


47.04 


30 


1.04 


2.263 


35.89 


50 


1.02 


2.363 


34.28 


100 


1.01 


2.440 


33.17 


oo 


1 


2.518 


32.13 
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Fig. 1. Density profiles of stellar polytrope for n < 5 (left) and n > 5(right). 
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Fig. 2. Trajectory of Emden solutions in (77, A)-plane. 
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Pc/Pe Pc/Pe Pc/Pe 

Fig. 3. Normalized specific heat per particle C^/N as a function of density con- 
trast Pel Pe near the critical polytrope indices n = 3(upper) and n = 5(middle), 
and large n cases (lower). Here, the normalized specific heat C* is defined by 

C v /(h 2 /GMr e )W 2 y( n - 3 / 2 l 
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Fig. 4. Stability /instability criterion for a system in contact with a thermal bath 
in (u, t>)-plane. The thick solid lines represent the trajectories of Emden solutions, 
while the thin-solid and dashed lines respectively denote the conditions (45) and 
(60). 
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Fig. 5. A schematic description of Carnot cycle. 
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